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Abstract 

The SPHC hydrodynamic code was used to evaluate the effects of orbital debris particle shape and orientation on 
penetration of a typical spacecraft dual-wall shield. Impacts were simulated at near-normal obliquity at 12 km/sec. 
Debris cloud characteristics and damage potential are compared with those from impacts by spherical projectiles. 
Results of these simulations indicate the uncertainties in the predicted ballistic limits due to modeling uncertainty 
and to uncertainty in the impactor orientation. 

In support of the Columbia Accident Investigation Board (CA1B), the NASA Johnson Space Center (JSC) Safety 
and Mission Assurance (S&MA) Directorate contracted with the Institute for Defense Analyses (IDA) to conduct a 
systematic review of potential causes of failure, including impacts by meteoroids or orbital debris (M/OD) (Ref. 1). 
One of the central findings in the IDA review was that NASA's critical M/OD risk predictions contain a number of 
significant input uncertainties. One of the largest of these uncertainties appears due to the lack of non-spherical 
shape considerations in NASA’s orbital debris environment and penetration models. As a first step in correcting this 
lapse, the report suggested that NASA perform a sensitivity analysis for the expected range of effects on damage 
considering spherical vs. non-spherical impactors. 

The latest version of the NASA orbital debris model is ORDEM2000, which was released for use in May 2002. 
This model utilizes updated in-situ impact data and ground-based radar data to form an empirical model of the 
current orbital debris flux, with predictions for flux growth in the out-years based on the NASA EVOLVE model. 
The Satellite Breakup Model (SBM) is the component of EVOLVE that derives an area-lo-size relationship for 
small orbital debris particles (Ref. 2) based on radar cross section (RCS) measurements from the Haystack and 
Goldstone stations. The debris particle “size” distribution is stated in terms of characteristic length, Lc. The SBM 
area-to-size relationship assumes that particles with Lc below i.66 mm are cubes, whereas particles above this size 
become increasingly “potato chip-” or “flake-” shaped, e.g., a particle with an Lc of 5.3 mm has a length-lo- 
thickness ratio of 3. These shapes are consistent with the RCS-lo-size conversion from the Haystack radar data and 
the mass-to-size relationships from the Spacecraft Orbital Debris Characterization Impact Test (SOC1T) fragment 
measurements (though not all SOCIT fragments had this type of shape). 

In order to examine the sensitivity of impact damage to both the shape of an orbital debris particle and its orientation 
at impact, we made use of the SPHC hydrodynamic code. This code implements the smooth particle hydrodynamics 
method in simulating impacts, and was previously compared to other codes and to the predictions of several familiar 
penetration equations (Ref. 3). For this study the code was run on Pentium IV ™ desktop computers under the 
Windows 2000 ™ operating system. The runs were fully three-dimensional simulations, reserving memory for up 
to 210,000 SPH particles. The projectiles and bumper layer were 2024 aluminum, and the backwall was 2219. The 
Whipple shield configuration was a 1.6-mm bumper, 12-cm standoff, and a 3.2-mm backwall, the same as that used 
in Ref. 3. 

We created an algorithm to assign length and thickness proportions to our impactors, based on data from SBM, for 
input values of Lc. We selected five impactor orientations to examine the effect of this parameter on the resulting 
debris clouds and damage to the backwalls. Sketches of our impact orientations are shown in Figure 1. The impact 
speed in all simulations was 12 km/s, and the velocity obliquity was 1.6 deg off the normal - this was done in order 
to break any gridding symmetries. 
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Figure 2 shows the simulation setups for a sphere and a flake, both with Lc = 0.6 cm. The sphere’s mass was 0.314 
g, and the Hake’s was 0.1 15 g. The Hake setup corresponds to the “Face A-li (45-45)” example in Figure 1 . Figure 
3 compares the debris clouds of these simulations at 8 microseconds into the events. The Hake's debris cloud lacks 
the low-density lead element of the sphere’s, and has a concentration of dense material in an arc at approximately 
the vertical locaiion ol the sphere’s central fragment mass. 1 Ins arc is aligned along what were the normals to the 
flake’s largest (square) faces. Nestled in the center of the arc is a barely-fractured remnant of the trailing corner of 
the Hake. The development of this arc structure along the Hake face normals is seen in debris clouds from other 
orientations as well. Figure 4 compares the backwall damage due to these projectiles at the 100 microsecond point. 
While both projectiles fail the backwall, the damage characteristics are distinctly different. The sphere produces a 
collection of smaii perforations near the center of the backwall, several of whose margins have coalesced iu make 
what should develop into a set of jagged petals at later times. The Hake produces a lenticular rip in the backwall, 
aligned with the dense debris arc and flake face normal line. Such a linear hole might cause the backwall to be more 
subject to an "unzipping” failure at later times if it happened to fall near the wall’s principle stress axis. At this 
stage the flake’s hole size is several times that of the sphere’s for this orientation, but whether this difference is 
maintained at later times is unknown. 


Hu and Schonberg (Ref. 4) reported that non-spherical impactors were much more damaging than spherical 
impactors of the same mass. We concur in this finding, since the case above shows that a square flake can produce 
greater damage than a spherical impactor of over 2.5 times the mass. However, the quantity observed in the debris 
population is not mass, but RCS, which is directly related to Lc. We prefer to develop ballistic limits based on this 
observed variable, but owing to the constraints of time, we can only report results obtained at 12 km/s, comparing 
flakes to spheres. 


Accordingly, we varied Lc for each of the five orientations shown in Figure 1, and recorded the Lc values for which 
the projectile produced four classes of backwall damage: clearly visible through-holes; spallation of SPH particles 
with acompanying backwall fracture, but no discernable through-holes; fractured back surface particles, but no 
spallation; and plastic deformation without fracture of the back surface. The Lc values, corresponding masses, and 
resulting damage are shown in Table 1. The Lc values in this table define a ballistic limit band applicable to these 
non-spherical particles at an impact speed of 12 km/s. The smallest values of Lc resulting in clearly v isible 
perforations define an upper bound to this band, each applicable to its specific orientation. The largest values of Lc 
resulting in plastic-deformation-only define the lower bound of the band, again for each sped flic orientation. For 
each orientation the shift from fracture-only to production of through-holes constitutes a transition depending on the 
specifics of the simulation model and code -a modeling uncertainty. Considering all the orientation cases together, 
we observe a spread in the Lc values for the various types of damage that is greater than the spread for any given 
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modeling uncertainty. By defining this wider band, we indicate an inherent uncertainty attached to any ballistic 
limit prediction that makes use of the SBM shape model. Thus, this uncertainty contributes to the damage 
uncertainty attachable to EVOLVE, and thence to the uncertainty in ORDEM2QOO. 


Future work will extend the velocity range covered down to at least 7 km/s, and up to 15 km/s. Statistical methods 
will be applied to refine our statements about the uncertainty bands, and any variation of these bands with velocity. 
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Table 1: Flake impact aspect run matrix. Velocity 1 2 km/s, obliuuitv 1.6 de». 


“Edge-on” Cases 

Lc 

Mass 


(mm) 

(g) 

Rx = 0 Ry = 0 Rz = 0 

4.0 

0.039 


4.5 



5.0 

0.071 


5.5 



6.0 

0.115 

Rz = 45 

4.0 

0.039 


4.5 



5.0 

0.071 


5.5 



6.0 

0.115 

Rz = 90 

4.0 

0.039 


4.5 



5.0 

0.071 


5.5 



6.0 

0.115 

“Corner-on” Cases 



Rx = 0 Ry = 45 Rz = 45 

4.0 

0.039 


4.5 



5.0 

0.071 


5.5 



6.0 

0.115 

Rz = 90 

4.0 

0.039 


4.5 



5.0 

0.071 


5.5 



6.0 

0.115 

Sphere Cases 

4.0 

0.093 


4.5 

5.0 

5.5 

6.0 


0.314 


Plastic Frac Spall Thru 


X 


X 


X 


X 


X 


X 


X 


X 
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Figure 1: Impact flake geometry and impact orientations used in the simulations. 



Figure 2: Sphere and flake setups, 3-D views. 
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Figure 3: 2-D views of debris clouds of sphere (left) and flake (right). 
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Figure 4: 3-D views of backwall damage by the sphere (left) and the flake (right); view is lrom 
behind the backwall. Dark blue material is fractured; light blue is plasticly yielded; green is 
undisturbed solid. Through-holes are visible as white. 
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